Comparison of averages of flows and maps 
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It is shown that in transient chaos there is no direct relation between averages in a continuos 
time dynamical system (flow) and averages using the analogous discrete system defined by the 
corresponding Poincare map. In contrast to permanent chaos, results obtained from the Poincare 
map can even be qualitatively incorrect. The reason is that the return time between intersections 
on the Poincare surface becomes relevant. However, after introducing a true-time Poincare map, 
quantities known from the usual Poincare map, such as conditionally invariant measure and natural 
measure, can be generalized to this case. Escape rates and averages, e.g. Liapunov exponents and 
drifts can be determined correctly using these novel measures. Significant differences become evident 
when we compare with results obtained from the usual Poincare map. 



Extensive investigations of chaotic systems in recent 
years have demonstrated the great importance of tran- 
sient chaos, due mainly to its connection with transport 
phenomena and chaotic advection [Q, possibly as- 
sociated with chemical reactions || . In most chaotic sys- 
tems it is sufficient to know the intersection points of 
the trajectories with a chosen surface P, the so-called 
Poincare surface. In case of iV-dimensional phase space 
P is N — 1-dimcnsional. Using a coordinate system on 
P, and finding the connection between the successive in- 
tersections x„ and x n+ i the Poincare map (PM) can be 
constructed as 



f(x n ). 



(1) 



The behavior of the system can then be studied by iter- 
ation of this map. The advantages of the use of PM are 
(i) it is discrete, (ii) it has smaller dimension. Its disad- 
vantage is the absence of the close connection between 
the number of intersections n and the time t, since the 
return time r between two intersections depends generi- 
cally on where a trajectory intersects. One can keep this 
information by completing the PM with the equation 



t n +i = t n + r(x„) 



(2) 



We call this extended map the true-time Poincare map 
(TPM). 

Usually, one reduces to the PM by the following ar- 
gument: The total time after n iterations is given by the 
sum of the corresponding return times t(x). It is gener- 
ally assumed that for large n and for typical trajectories 
the terms in the sum can be replaced by their average 
over the invariant density pp of the map. The sum then 
becomes a product H 



t = n(r), (t) 



dccpp(x)r(x), 



(3) 



Based on this connection, averages of the map (using n 
for time) and the flow (using the real time t) would be 
simply related by a time scale. This is explicitly shown 



for general averages in case of permanent (non-transient) 
chaos [Q. 

We demonstrate in this paper that, in contrast to per- 
manent chaos, the situation is quite different for transient 
chaos. Not only (r) in Eq. (||) should be changed, but 
averages of the map and the flow (or of the TPM rep- 
resenting it) are not any more related by a time scale. 
The situation is somewhat reminiscent of the case when, 
instead of simple averages, the decay rates of correlations 
are considered. Even in permanent chaos, these show a 
discrepancy in non- ideal situations To proceed in a 
correct manner we must start with the TPM, which con- 
tains all the information needed for the longtime behavior 
of the system and from which we can derive all necessary 
formulas. Finally we compare these with the correspond- 
ing ones of the PM by setting r(x) = (r). The use of 
the PM is sufficient if the results do not change. 

It is convenient to initiate the trajectories by inserting 
particles on P with an input current density pi n (x,t). 
Let pi n (x,t) = for t < 0. Since a trajectory leaving 
P has either been initiated there or has intersected P 
previously, we obtain for the normal component pp(x,t) 
of the current density on P 



pp = (Cp P ) + pi-, 



(4) 



Here C is the Frobenius-Perron operator of the TPM 
which is defined by 



(£g)(x,t) = J dx'5(x 



f(x')). 9 (x',t-r(x')). (5) 



Quite often the motion in one direction — the unstable 
one — depends at most weakly on the others. Choosing 
a coordinate system in which x is taken along this di- 
rection evolution of x can be well approximated by a 
one-dimensional map x n +\ = f(x n ). (The price paid is 
the non-uniqueness of / _1 ). This happens, for example, 
in strongly dissipative systems and in those analogous to 
Baker type maps. In such a situation Eq. (Q) remains 
valid if x is replaced by x and f by /, and projecting the 
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densities onto the unstable direction. For simplicity we 
restrict our attention to this one-dimensional case. 

First we compute the quasistationary distribution. We 
assume that the system will become quasistationary after 
some time, i.e. that the distribution decays exponentially 
but all relative weights remain constant. Normalization 
of the distribution leads to the time-independent con- 
ditionally invariant density pp{x). We make the 
ansatz pp(x,t) = pp{x)e~ Kt , where n is the escape rate, 
and obtain from Eq. (^) the selfconsistent equation for 
the conditionally invariant density 

p P (x) = Jdx'5(x - f{x'))e KT ^ PP {x') . (6) 

Here / is the range of the values of x. 

For example, let / be the tent map with a possible 
opening: f(x) = x/clq if x < ao, fix) = (1 — x)/a\ if 
x > 1 — a%, where ao + &i < 1. Furthermore let t[x) 
be piecewise constant: t(x) = tq if x < ao, t{x) = t\ 
if x > 1 — a\. Escape occurs for ao + a\ < 1, when 
the trajectory leaves the Poincare surface in the interval 
x G (ao, 1 — ai). The smooth, non- negative solution of 
(||) is now pp(x) = 1. By chance pp(x) does not depend 
on r(x), but the equation for n does: 



a e KT ° + aie h 
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and in nonlinear maps pp{x) also depends on it. 

For a general treatment we write the formal solution 
of Eq. (H) as 



pp(x,t) = [(1 - C) 1 p in ](x,t) , 



(8) 



We continue with Laplace transformations in time, since 
the generalized operator can be written as (Cg)(t) = 
f *dt' L(t-t')g(t'), with L(At)g{x,t') = Jjdx'6(f(x') - 
x)S(t(x') — At)g(x' , i'). Its Laplace transform (denoted 
by ~ ) is Cg = Lg and Eq. (ji|) yields 



pp = (1 - L) 1 p in 



(9) 



Considering s as a parameter we can use the eigenfunc- 
tions satisfying 



L(s)tp m (s) = X m (s)ip m (s) 



(10) 



to expand p in as p in {s) = ^o° a m (s)ip m {s). Inverse 
Laplace transformation gives 



e+ioo 



Pp(x, t) = -^7 / ds V - — " 

27TI J c - ioo n 1 ~ *™(S) 



a m (s)cp m (x, s) . 



Each value of s for which A m (s) = 1 with some m gives 
a pole in the integrand and a term e st in pp(t). There- 
fore, the leading asymptotic time dependence is e _Kt , and 
hence the escape rate k is determined by the position of 
the leading pole, i.e. 



k = —So, so = max{s with A m (s) = 1 and s real}, (11) 

rn 

where we assume for simplicity that s is maximal for 
m = 0. (s must be real, otherwise pp{t) could not re- 
main positive for all t.) Eq. ( |To| ) together with Eq. ( pd| ) 
corresponds to (^|), however, here we have obtained the 
result and the decay 



pp(x,t) 



-e Kt a tp (x, 



(12) 



for large times without prescription of e~ Kt . 

Typical long time averages are determined by trajecto- 
ries in the following manner. The quantities in question 
consist for each trajectory of a sum of terms A(x) that 
are added whenever the trajectory intersects P, i.e. we 
need averages of expressions of the form A(xi) with 
xi = f l (x ). We consider some examples: If we are inter- 
ested in the number of intersections n, we set A(x) = 1, 
for the total time we set A(x) — r(x). If we are inter- 
ested in the leading Liapunov exponent describing the 
exponential deviation of infinitesimally close trajectories 
we need the logarithm of the full derivative of /, i.e. we 
must set A{x) = In | f'(x) |. 

We now outline the case of the ordinary PM . The 
average of ^2[A(xi) is obtained by taking into account 
the contribution of all trajectories present after n it- 
erations and dividing this by the weight of these still 
present trajectories. Thus we arrive at Jj dx Jj dxo S(x n — 
a: )I]/!=o 1 A(xi)p in (x Q )/(f I dx p P (x,n)). The denomina- 
tor is an integral of the density pp(x, n) = L dxo S(x n — 
x))pi n (xo). Correspondingly we can write the numer- 
ator as an integral of a 'weight density' 071(3;, n) — 

Jj dx S(x n - X) ^27=0 A(xi)p in (xo). 

The TPM requires two modifications: (i) instead of n 
we take a time interval of length i, (ii) an additional sum 
over the number of intersections n is necessary, since n is 
generically different for different trajectories. This yields 
for the weight density 



<TA(x,t) = / dx ° / dto8(x n - x) 

n=0 ^ J 



(13) 



n—1 n—l 

x S(t - t - t{x 3 )) X A(xi)p in (x ,to) , 

j=0 1=0 

where S(t — to — ^2" = q T ( x j)) demands that the trajec- 
tories intersect n times during time t — to. The average 
value of A(xi) is given asymptotically as 

t(A) = t lim 1 IidxaAip^) ^ 
t'—>oo t' Jj dxpp(x,t') 

The index (i signalizes that asymptotically this average 
does not depend on pi n but on the natural measure p of 
the TPM. This will now be explained. 
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To compute a a we note first that 
a a = C{a A + Ap P ) . 



(15) 



TPM differ significantly, (although their dimension Do 
is the same.) 

The Liapunov exponent can be written as 



This is analogous to Eq. @. Defining T = (1 — C) 1 we 
obtain 



a A — TCATp 



(16) 



We introduce T by (Tg)(t) = J*dt'T(t - t')g(t'), sim- 
ilarly to the connection of C and L. We can write the 
Laplace transform of ( fhf| ) in terms of the adjoints L + 
of L and T + of T. For the backward transformation of 
this expression we need the eigenf unctions of L + . The 
solutions iprn(s) of L + (s)ip m (s) = AJ^(s)^> m (s) are func- 
tionals p^jl^] (and can be approximated with strongly 
oscillating functions), due to the fractal nature of the in- 
variant set. We insert an expansion 1 = X)o° ^m^m and 
observe that, for large t, the most important terms occur 
when poles induced by T + *(s) and T(s) coincide. Thus 
we obtain for large t 

/ dx cr A (x, t) = -—J! — dx ipo{x, -n)A(x)<po(x, -«). (17) 

J I A o \ K ) J I 

(Note that both i/jq(—k) and lpo(~k) are real.) The pref- 
actor of Eq. ([I?]) can be expressed by setting A(x) — t(x) 
because of the relation, valid for large times, (r) M = 1. 
Asymptotically we obtain 



{A) » = V)L 



dx ipo(x, — K)ipo(x, — k)A(x) , (18) 



where (r) = U dxipo(x, — K)(po(x, — k)t(x) with the nor- 
malization Jj dxipo(x, —K,)(po(x, —k) = 1. Since these re- 
lations are valid for every observable A, we identify the 
natural measure ||13| for infinitesimal intervals as 




FIG. 1. Fractal distributions of the natural measures when 
modeling f(x) by the open tent map (ao = ai = 0.475) left: 
normal Poincare map (scale n), To = Ti = 1. right: true-time 
Poincare map (scale t), to = 1, ri = 0.1. 

fi([x, x + dx)) = —■ip (x, —K)ipo(x, —k) dx . (19) 
v / 

A comparison between the natural measure of the PM 
and the TPM for the tent map (Fig.l) shows obvious 
differences. It is clear that the repeller of the PM and 



A 



Liap 



dpLp In | f'{x) 



(20) 



For A(x) = 1 we find n(t) = i(l) M = t/(r), thereby 
< = (r)n . 



(21) 



If we compare Eq. (yj), Eq. (|21j) and Eq. fl20J) with Eq. 
(^) we see a complete analogy, except that r has to be 
averaged over a quantity of the TPM. This means that 
the flow of a repelling system is well represented by its 
TPM and not by its PM. 

The leading Liapunov exponent for the repelling tent 
map is 

1 a expKr ln(a " 1 ) + oi expKTi ^(a^ 1 ) 

^Liap = 73" ^~~ZZTZZ — I I -_ • 



ao exp ktq + ai exp kti 



In this example we see the irrelevance of r for n = 
(no transient chaos), so that (r) only sets the time scale. 
On the other hand, for n > ALi ap is not invariant to 
changes of To relative to t%, proving again that r is a 
relevant quantity in transient chaos. 

The comparison of the behavior of the flow and the 
map shows that the occurrence of criticality can change 
when turning from the map to the real system. A state 
of a system is called critical if the natural measure is 
concentrated on a subset of the repeller, while the in- 
variant measure is distributed on the whole repeller. In 
such situations there are two conditionally invariant mea- 
sures with different escape rates |l5|,[l6| , p^| . Here we 
use the piecewise parabolic ID map|15[ that is defined 
on the interval [0,1] by its inverse branches f^ 1 (x) — 
(x + d-x(l-x))/2R (lower branch), f^ix) = 1 - 
ff 1 {x) (upper branch) and choose t(x) = 1+t x -(x— 1/2). 




FIG. 2. The leading escape rate of the piecewise parabolic 
map as a function of the derivative of t(x),t x . 

When increasing t x at a certain value the escape rates 
change order of the modulus. This is evidently a break- 
point in Fig. 2, which shows the leading (smaller) one of 
the escape rates. Above that point criticality disappears. 
This will be explored in more detail elsewhere. 
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If P has a periodic structure it can be reduced to a 
unit cell with periodic boundary conditions. In this case 
both PM and TPM can be separated into a reduced map 
(which maps the unit cell into itself) plus a shift A(x) 
describing the transit between the cells 0,^8) ■ Such sys- 
tems can be characterized by the drift speed and diffusion 
coefficient. If particles can be lost from the point of view 
of diffusion by absorption, chemical reaction or escape in 
directions transverse to the extension of the system, we 
refer to transient diffusion p9|-pl|. We set A = A and 
we obtain a shift density cta in analogy to the procedure 
above: a a — T£ATp in . We then determine the drift 
speed as the normalized shift per time 

average of (Sft) } 1 , . , ,„„, 
V =hm % = 5)< A >"- (23) 

As an example of the essential role of averages over 
the natural measure we consider a diffusive system on 
a one-dimensional lattice and assume that the reduced 
map is the tent map. We assume furthermore a micro- 
scopic process determining whether and by how much a 
particle jumps to the left or to the right on the lattice. 
We consider a lattice with period 1 and a change of the 
coordinate by ±1 depending on the location in the subin- 
terval [0, 1], namely A(x) = — 1 if x € [0, a ], A(x) = 1 if 
x 6 [ai,l]. To calculate the average speed or the diffusion 
coefficient an average over long trajectories is required. 
The average speed is 

I ai exp kti — ao exp ktq , _ , , 

v = —— . (24) 

(t) eti exp kt\ + a,Q exp kto 

Again if k = (non transient chaos) the return time r 
sets the time scale only, and PM and TPM give the 
same result. But for k > (transient chaos) even the 
sign of the speed can change when computing it with the 
usual PM, i.e. when setting t(x) = (r) = const. Results 
for the diffusion coefficient will be published elsewhere. 

In conclusion we have shown that the return time r 
i.e. the time between two successive intersections on the 
Poincare surface P, is a relevant quantity in transient 
chaos. The usual Poincare map does not reflect the long 
time averages of the flow satisfactorily and can even be 
completely misleading. The solution is to use a true- 
time Poincare map TPM and its generalized Frobenius- 
Perron operator (^|), where we can also define condi- 
tionally invariant measure and natural measure. Escape 
rate, Liapunov exponents, drift speed, etc. depend signif- 
icantly on t(x) and are described correctly only by using 
the TPM. Therefore the necessary generalization of the 
normal Poincare map is the true-time Poincare map if 
the system in question is a repeller. 
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